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TECHNICAL NOTE 2156 


THEORETICAL CALCULATIONS OF THE LATERAL FORCE AND YAWING 
MOMENT DUE TO ROLLING AT SUPERSONIC SPEEDS FOR SWEPTBACK 
TAPERED WINGS WITH STREAMWISE TIPS ' # 

SUPERSONIC LEADING EDGES 
By Sidney M. Harmon and John C. Martin 

SUMMARY 


Theoretical results are obtained, by means of the linearized theory, 
for the lateral force and yawing moment due to combined angle of attack 
and roll. The thin wings considered are tapered and have leading anH 
trailing edges which are each swept at a constant angle; the tips are 
parallel to the center section. Further restrictions are that the leading 
edges are supersonic (the velocity component normal to the lea ding edge is 
supersonic) and the Mach lines from either tip do not intersect the remote 
tip. The trailing edges may be either supersonic or subsonic provided 
that the taper ratio is equal to or less than one. 

Design curves are presented which permit rapid estimation of the two 
stability derivatives, rate of change of lateral force and yawing moment 
with roll Cy p and C np , respectively, for given values of aspect ratio, 

taper ratio, Mach number, and leading-edge sweep. These design curves 
are given relative to the body axes with the origin located at the pro- 
jection of the leading edge of the tip section on the center section. 

The conversion to stability axes with moments about an arbitrary, axis can 
be accomplished by use of a simple formula given herein. 


INTRODUCTION 


One of the important derivatives required in the prediction of the 
lateral stability of airplanes is the rate of change of yawing moment 
with roll. Theoretical results for this derivative Cnp and also for 

the derivative expressing the rate of change of lateral force with 
roll CYp at supersonic speeds have been reported in a number of papers 
for thin wings having particular plan forms. These papers include 
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references 1 and 2 for triangular plan forms, reference 3 for rectangular 
wings, reference h for wings of zero taper ratio having sweptback, sub- 
sonic leading edges, and reference 5> for sweptback tapered wings with 
streamwise tips and subsonic leading edges (the component of velocity 
normal to the leading edge is subsonic). 

In the present paper, the linearized theory is used to obtain the 
derivatives C n p and Cjp for thin wings at supersonic speeds. The 

wings considered have conventional taper ratios, the leading and trailing 
edges are each swept at a constant angle, and the tips are parallel to 
the center section. The results are valid only for a range of super- 
sonic speeds for which the leading edge is supersonic (the component of 
velocity normal to the edge is supersonic). The trailing edges may be 
either supersonic or subsonic provided that subsonic edges dcf not affect 
the tip edges. A further restriction is that the Mach lines from either 
tip do not intersect the remote tip. 
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SYMBOLS 

9 

rectangular coordinates (figs. 1 and 2(a)) 

induced flow velocity along Y-axis 

undisturbed flight velocity 

angular velocity about the X-axis 

density of air 

speed of sound 

stream Mach number (V/a) 


Mach angle 



wing angle of attack in steady flight, radians 
sweepback angle of wing leading edge 


m' = B cot A 
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s n 
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A' - BA 


disturbance-velocity potential on upper surface 
of airfoil 

ratio of tip chord to root chord (taper ratio) 

wing span 

total wing area 

distance to tip edge 

aspect ratio 


x C g distance from wing apex to center of gravity in 

root chords measured positively as indicated 
in figure 2(b) 

F induced suction force on wing tip per unit length 

of tip 

G function of x given by equation (l) 

Cj 8 Rolling moment 
| V 2 Sb 


G y = 


Lateral force 
| V 2 S 


C n = 


Yawing moment 

| V2sb 


% = | 

( Gn p)l 



with respect to body axes with origin located 
x cg root chords from leading edge of center 
section 


with respect to body axes with origin located 
at projection of leading edge of tip section on 
center section (fig. 2(a)) 
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Subscripts: 

TE 


with respect to stability axes with origin located 
Xgg root chords from, leading edge of center 

section 


with respect to body axes 


refers to wing trailing edge 


SP 


contributions caused by angle of attack and 
rolling motions, respectively 


ANALYSIS 

Scope 


The wing configurations considered in the analysis are shown in 
figure 1. These wings have streamwise tips and sweptback leading edges, 
although the trailing edges may be either sweptback or sweptforward. 

The results are valid only for a range of supersonic speeds for which 
the leading edge is supersonic. The trailing edge may be either sub- 
sonic or supersonic provided that the taper ratio is equal to or less 
than one. For subsonic trailing edges the Kutta condition is assumed 
to be satisfied along the trailing edges. A further restriction on 
these wings is that the Mach lines from either tip do not intersect the 
remote tip. The analysis is based upon linearized theory; therefore, 
the results are subject to the restrictions of this theory. 

The derivation of the stability derivatives is made initially with 
reference to the body axes which are fixed in the wing with the origin 
located at the projection of the leading edge of the tip on the center 
section (fig. 2(a)). The transformation of the stability derivatives 
to a system of stability axes (fig. 2(b)) is discussed in the section 
entitled "Derivation of Formulas." 
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Basic Considerations 

The stability derivatives are determined ‘from, integrations of the 
farces and moments over the wing. If skin friction is neglected, these 
forces and moments arise frcm. pressures on the wing surfaces and from 
suction forces along the edges of the wing. A lifting wing with stream- 
wise tips or subsonic leading edges will induce suction forces along the 
subsonic leading edges and along the streamwise tips. These forces arise 
because of the flow from the bottom to the top surfaces. This analy sis 
is concerned only with the derivatives Cy and C n . For these deriva- 

r r 1 

tives with respect to the body axes, the lateral forces and yawing moments 
are produced only by suction forces along the edges of the w ing . Sinc e 
in this analysis the leading edges are supersonic, no suction forces are 
induced along these edges j thus, the determination of the derivatives Cy-p 

and Cjjp in the cases considered has involved only the evaluation of 

unbalanced suction forces along the wing tips. The Kutta condition is 
assumed to be satisfied along the subsonic tr ailing edges. 

The suction forces along an edge of infinitesimally small thickness 
can be evaluated theoretically by a method suggested in reference 6 for 
incompressible flow and modified for compressibility effects in refer- 
ence 7 . If this method is used, the suction force varies as the square 
of the induced flow velocity around the edge. In general, in the linear- 
ized theory the induced velocity varies directly as the magnitude of a 
given motion. Thus for constant angle of attack, the induced velocity 
around the edge is proportional to aj and in a rolling motion, the 
induced velocity .is proportional to p. For constant angle of attack 
and constant rate of roll, the suction forces vary as the squares of a 
and p, respectively. These forces are in the opposite direction and 
of the same magnitude j thus no lateral force or yawing moment is devel- 
oped. For combined angle of attack and roll, unbalanced suction forces 
occur that vary as the cross product of a and p; consequently, these 
forces or the associated stability derivatives cannot be obt aine d by 
superposing separately the components due to each motion. The present 
analysis takes into account this general condition of simult ane ous 
ro llin g and vertical motions. 


Derivation of Formulas 

According to the methods of references 6 and 7, the induced surface 
velocity normal to the wing tip is expressed as 


C v n) s n — > 0 


G 


( 1 ) 
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where G is a function of x only and s n is the distance to the tip 
edge. Then the suction force per unit length of tip is 

F = npG 2 (2) 


for streamyd.se tips. 

The induced surface velocity normal to the tips of a wing simul- 
taneously in roll and at a given angle of attack a is given by 


( v n) 


a+p 



(3) 


where the quantities are evaluated in the immediate vicinity of the wing 
tip, that is for s n — >0. The potential functions 0 a and 0p for use 

in equation (3) were obtained from table I (regions 3 and 5) of refer- 
ence 8. By utilizing these expressions for the potential functions and 
by substituting for v n in equation (l), the function G is determined. 
The suction force per unit length along the wing tip F is then obtained 
from equation (2). 

If the origin of the body axes is taken at the projection of the 
leading edge of the tip on the center section (fig. 2(a)), the non- 
dimensional forms of the expressions for the derivatives are as follows: 



k i F dx 

J Tip 

pVSbp 


(U) 



k F Fx dx 

J Tip 

pVSb2p 


(S) 


where the integrals are evaluated along both wing tips. 

Formulas for the derivatives Cy^ and are given in the 

appendix. These formulas are for the two different plan-form configura- 
tions which occur, that is, for the cases where the foremost Mach line 
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the wing trailing edge or the tip. (See 

is readily transformed to a system of body 

axes with its origin at an arbitrary point along the center section. 
Thus 3 if the origin of the axes is located XQg root chords from the 

leading edge of the center section, the following equation results 


from the center section cuts 
fig. 3.) 


The derivative 


(NX 




A(1 + X) 



( 6 ) 


The nondimensional stability derivatives Cy^ and which 

were derived with reference to body axes with the origin located at 
the projection of the leading edge of* the tip on the center section 
(fig. 2(a)), may be transformed to apply to stability axes with the 
origin at an arbitrary point along the center s.ection. (See fig. 2(b).) 


The transformation formulas for conversion from body axes to sta- 
bility axes are given in reference 9 and also in table II of reference 3 
where terms of the order of a£ are omitted. A consideration of these 
formulas and the derivatives involved therein shows that, to the first 
order in a, the conversion formulas from body axes to' stability axes 
are 



(7) 



1 


Gn p - ad? = (Cnp)! + 


2x cg 
A(1 + X) 



T P 



( 8 ) 


where the primes refer to the stability axes with the origin located 
along the center section at an arbitrary distance x cg root chords from 

the leading edge. 

The derivative which appears in equation (8), if obtained to 

the first order in a, has the same value either with respect to body 
axes or with respect to stability axes. Values for Cjp for use in 

equation (8) may be obtained from reference 8 for all the wing configu- 
rations considered in this paper, except for a special case. This 
special case occurs if the Mach line from the leading edge of a tip 
section intersects the remote half -wing coincidentally as the Mach line 
from the leading edge of the center section intersects the tip section. 
(See fig. 1(b).) 
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RESULTS AND DISCUSSION 


The results of the calculations are presented in figures 4 and 'j in 
the form Cyp^a and ^Cnp^/Ba plotted against B cot A for various 

values of BA' and X. The derivatives Cyp and (Cnp)y were computed 
from equations (Al) and (Af>) and from equations (A2) and (A6), respectively. 


If in figures U and 5 the Mach number and B are constant, the data 
in these figures indicate directly the variations, respectively, of the 
derivatives Cyp and with the angle of sweep of the leading 

edge, aspect ratio, and taper ratio. With B held constant, the data in 
figures U and $' indicate that as the leading-edge sweep angle is increased 
(decreasing values of {3 cot A) Cyp increases positively and (j^np)y 

increases negatively. The data in figures U and 5 also indicate that, 
with B, X, and A constant, Cy p and (C np j 1 increase in magnitude 

with decreasing aspect ratio and, with B, A, and A constant, Cyp 


and 



increase in magnitude with increasing taper ratio. 


Figure 6 presents some illustrative results for the derivative C np ' 

with respect to the stability axes. These data illustrate the variations 
in C n p' with Mach number, aspect ratio, sweepback, and taper ratio. 

For these illustrative results the origin of the stability axes is 
assumed to be at the leading edge of the center section (x^ = 0). 

Since the origin is at the leading edge of the center section, the sweep 
effect shown is considerably greater than would probably be encountered 
in an actual airplane design where the longitudinal location of the 
center of gravity would probably be near the wing center of pressure. 

The effect on C n p‘ of shifting the origin of the stability axes is not 

shown in figure 6 but is readily obtained from an examination of equa- 
tion (8). From equation (8) shifting the origin of the stability axes 
rearward (increasing Xgg^ is noted to result in a change in the yawing- , 

moment coefficient equal to 


2xc g (GYp) 

A(1 + X) 

Thus, for the usual case in which Cy p is positive, an increase in x cg 
contributes a positive increment to C^. The data for C n ' given in 

r * 

figure 6 were computed from equation (8), in which the values for Cyp 
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and ( c npj-]_ were obtained from figures k and 9, respectively, and the 
values for Cj were obtained from figure 10 of reference 8. 

The results of the present investigation have been derived on the 
assumptions of zero thickness and small disturbances. Theoretically, 
infinite lateral perturbation velocities are induced along the tips. 

These velocities cause infinite perturbation pressures which, when acting 
on an edge of zero thickness, yield a finite suction force. In the 
actual flow, however, the infinite pressures are never realized; but the 
drop in pressure tends to be compensated by the increased force which 
results from the finite thickness of the actual side edge (wing tip). 
These points are discussed in slightly more detail in reference 2. 


CONCLUDING REMARKS 


General expressions have been obtained by means of the linearized 
theory for the lateral force and yawing moment due to rolling for swept- 
back tapered wings with streamwise tips and supersonic leading edges. 

The trailing edges are swept and may be either supersonic or subsonic 
provided that the taper ratio is equal to or less than one. A further 
restriction on these wings is that the Mach lines from either tip do not 
intersect the remote tip. The analysis is based upon linearized theory; 
therefore, the results are subject to the restrictions of this theory. 

The results of the investigation were presented in the form of 
generalized design curves for rapid estimation of the two stability deri- 
vatives, the rate of change of lateral force and of yawing moment with 
roll CXp and Cn p , respectively. Some specific variations of the 

more important derivative C n p with Mach number, aspect ratio, leading- 
edge sweep angle, and taper ratio were also presented. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Ya., May 18, 19^0 
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APPENDIX 

FORMULAS FOR Cy p AND 


Formulas are given for Cy p and vr ^ 1 res P e ct to the body- 

axes with its origin located at the projection of the leading edge of 
the tip section on the center section for thin tapered wings with swept 
lea ding and trailing edges and with streamwise tips. The leading edges 
are supersonic and the Mach lines from either tip do not intersect the 
remote tip. These formulas for the two cases illustrated in figure 3 
are limited by the following conditions: 

B cot A, ^ 1 

B cot Ajrg = 1 or B cot Aipg < 0 


and 


BA > 


21 

1+1 


Case I - Foremost Mach line from center section cuts wing trailing 
edge.- For case I for which the foremost Mach line from the center sec- 
tion cuts the Wing trailing edge, the equations for CYp and ^Cn p ^ 


are 


°Xp = 


_ [9A'(1 + l)(m» - 1) - 81m'] 

9n(A'» ) 2 (m» - 1) 2 (1 + 1)3 


(Al) 


, n = -256aBl3m' [A» (1 + l)(m' - l) - lm«] 

' ^ 3n(A') 3 (m' - 1) 2 (1 + 1)^ 


(A2) 


When the leading edge is straight (m = °°), equations (Al) and (A2) 
reduce to the following expressions: 


gt. 


6UccX 2 [9A t (l + 1) - 81] 
9n(A') 2 (l + l) 3 


(A3) 
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fn \ 2S6oEX 3 [a*(1 + X) - x] 

lpjl 3n(A.) 3 (l + X)U 


(Alt) 


Case II — Foremost Mach line from center section cuts wing tip.- 

For case II for which the foremost Mach line from the center section 
cuts the wing tip, the equations for Cjp and ^Cnp^ are 




ka 


P 9"(A' ) 2 m« (m* + l) 2 


—(A* ) 3 (m* - l)(li{m« + 6) + 


l6m’X 


[3A» (3m' + 1) + 3A«m»X(3m» - l) _ 8m»X 2 

L 1 + x (1 + X) 2 (1 + \)3J 


(AS) 


liaB 


( Cl ‘ P ) 1 9n(A') 3 (m') 3 (m' t l) 2 

l6( m .) 2 A 2 [- 3(A ' )2(3m V 

L & * x ) 2 


(A* )^(m* - l) 2 (2m« + 1) + 


1) UA»m»X(3m» - 1) 12 (m* ) 2 X 2 

(1 + X) 3 (1 + \)U J 


(A 6) 


* When the leading edge is straight (m = ®), the expressions for 
Cjp and (pnp)]^ S iTen by equations (AS) and (A6), reduce to the same 

corresponding formulas as those obtained previously for case I. Thus, 
Cy p and ^Cnp^ for case II when the leading edge is straight are 

given by equations (A3) and (Al|), respectively. 
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(a) Sweptback leading and trailing (b) Sweptback leading edge and sweptforward 

edges. trailing edge. 


Figure 1.- Wing configurations considered in analysis. Supersonic leading 
edges j streamvise tips. Note that Mach lines from leading edge of 
center section may intersect either the tip or trailing edge and that 
Mach line from either tip does not intersect the remote tip. 


B cot A = 1, B cot Ajrg = 1 or B cot < 0, and BA > 

1 + ^ 
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Ob) Stability axes. (Corresponding body axes dashed for comparison.) 


Figure 2.- Systems of axes. 





(a) Case I. (t) Case II. 

Figure 3 • - Two types of plan-form configuration used for derivation of 
formulas for lateral-force and yawing-moment derivatives. 
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(a) Taper ratio \ = 0.25. 

Figure 5 .- Variation of ( c n p ) 1 | Ba ' wit]l sweepback-angle parameter for 

■various aspect-ratio parameters and taper ratios. Derivative is with 
respect to body axes with origin at projection of leading edge of tip 
on center section (fig. 2(a)). (See equation (8) for conversion to 
stability axes.) 
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(c) Taper ratio X = 0.75' 
Figure 5 .- Continued. 
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(a) Variation with Mach number. (b) Variation with aspect ratio. 

A = 53°; = 0.5. M = 2j A « 53°J >• = 0.5. 

Figure 6.~ Illustrative variation of C^ya with Mach number, aspect 

ratio, sweep angle, and taper ratio. Derivative is with respect to 
stability axes with origin at leading edge of center section. 
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